Becke-Johnson-type exchange potential for two-dimensional systems 
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We extend the Becke- Johnson approximation [J. Chem. Phys. 124, 221101 (2006)] of the ex- 
change potential to two dimensions. We prove and demonstrate that a direct extension of the 
underlying formalism may lead to divergent behavior of the potential. We derive a cure to the ap- 
proach by enforcing the gauge invariance and correct asymptotic behavior of the exchange potential. 
The procedure leads to an approximation which is shown, in various quasi-two-dimensional test sys- 
tems, to be very accurate in comparison with the exact exchange potential, and thus a considerable 
improvement over the commonly applied local-density approximation. 
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I. INTRODUCTION 

The advent of density-functional theory 1,2 (DFT) was 
followed by the development of approximate functionals 
for the exchange and correlation energy of many-electron 
systems. Significant advances were achieved beyond 
the local-density approximation (LDA) by generalized- 
gradient approximations, orbital functionals, and hy- 
brid functionals.- Those efforts focused almost solely on 
three-dimensional systems. However, already since the 
1970s systems of reduced dimensionality became of great 
relevance in solid state physics. The present nanode- 
vices consist of a large variety of low-dimensional sys- 
tems, where the many-body effects of interacting elec- 
trons need to be addressed. Of particular interest are 
two-dimensional (2D) structures including, e.g., semicon- 
ductor layers and surfaces, quantum Hall systems, and 
various types of quantum dots4 

Previous studies have shown that density functionals 
developed particularly for 3D fail when applied to quasi- 
2D systems^i By "quasi-2D" we mean physical sys- 
tems whose main building block is the quasi-2D electron 
gas. This class comprises most systems listed in the pre- 
vious paragraph and they can be reliably treated by a 
"pure" 2D approach, i.e., on a 2D grid. The physical 
reason for this is that the confinement potential along 
the growth direction (say z), is typically much stronger 
than in the two other spatial dimensions. This results in a 
quantum mechanical supression of the degrees of freedom 
along the z direction, and the system becomes effectively 
2D. Then, the influence of the surrounding host material 
is dealt within the effective-mass approximation mani- 
festing itself as an effective mass and a dielectric constant 
in the 2D Hamiltonian. 4 

Within the DFT approach, 2D many-electron prob- 
lems are most commonly treated using the 2D-LDA ex- 
change functional^ combined with the 2D-LDA correla- 
tion parametrized first by Tanatar and Ceperley^ and 
later, for the complete range of collinear spin polariza- 
tion, by Attaccalite and co-workers^ Despite the rela- 



tively good performance of LDA with respect to, e.g., 
quantum Monte Carlo calculations^ in semiconductor 
quantum dots, there is a clear lack of accurate 2D den- 
sity functionals to deal with diverse situations, espe- 
cially in the strongly correlated regime. Only recently, 
the transition from 3D to 2D has been analyzed,— and 
exchange-correlation energy functionals specially tailored 
for the 2D world have been P roposed.± 3 -±*»i5i±£»iL±S»i2 1 2a 
Also very recently, a DFT framework specially tailored 
for strongly correlated electrons has been developed and 
successfully tested for few-electron quantum dotsi2! 

In the present work, we aim at moving another step 
along the direction from 3D to 2D. We focus on approx- 
imations for the exchange potential, where our aim is 
to achieve an accuracy comparable to the exact result 
- available from the optimized-effective-potential (OEP) 
metho d 22 ' 23 or its approximations such as the Krieger- 
Lie-Iafrate (KLI) potential 2 ^ - yet simplifying the com- 
putational effort. We are able to find such an approxi- 
mation by (i) extending the framework of the well-known 
Becke- Johnson potential 2 ^ to 2D, and (ii) requiring the 
potential to satisfy the gauge invariance and correct 
asymptotic behavior. Our spin- and current-dependent 
exchange potential is indeed well comparable to the KLI 
and superior to the LDA in 2D quantum dots with vary- 
ing sizes and varying external magnetic fields. 



II. FORMALISM 



A. Spin-density functional theory 



Considering open-shell systems, (collinear) spin-DFT 
(SDFT) formalism 25 is one of the most useful. In this 
case, the total energy £ of a 2D system of interacting 
electrons is a functional of the two spin densities, p a (v) 
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with a =|, j: 



where 



+ / d 2 r v Q<J {r)p{v) + U[p]+ E xc \p hPl ] , (1) 
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Ts[ PhPl ] = E £/ ^ 
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is the non-interacting kinetic energy, the sum runs over 
occupied states, and N a is the number of electrons with 
spin a. Hartree atomic units (a.u.) are used, unless 
stated otherwise. Above, v 0a (r) is the spin-dependent 
external potential. The classical electrostatic (Hartree) 
interaction energy is given by 



d 2 r 



d 2 r 



, P(r)p(r') 



where 



p(r) = p T (r) +pi(r) 



(3) 



(4) 



is the total electronic density, which is the sum of the 
spin densities 



pA r ) = E I^M 



(5) 



with the sum running over occupied orbitals. Besides, 
E X c\P\ i Pi] is the exchange-correlation energy functional, 



which in practice needs to be approximated. The single- 
particle orbitals <pj a {r) in Eq. ^ are solutions of the 
Kohn-Sham (KS) equation^! 



— — + v sa {v) ) ip ja (r) = e ja ip ja {r) , (6) 



where j is a collective index for the one-electron quantum 
numbers, except for the spin. The effective single particle 
potential for spin a is given by 



v sa (r) = V 0(T (r) + u H (r) + v xca (r) 



with the Hartree potential 



v H (r) = J d 2 r 



, P(r') 
|r — r' 



(7) 



(8) 



and the exchange-correlation potential 

5E xc [p h pi\ 



5pa(r) 



(9) 



B. Optimized-effective-potential method 

In the following, we restrict ourselves to orbital func- 
tionals. A standard example of this type of functionals 
is the exchange energy functional of the form 



Ex [PhPi] 



EE 
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(10) 
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which resembles the Fock term of Hartree-Fock theory, 
here evaluated with KS orbitals. This is the so-called 
exact- exchange (EXX) energy functional. It should be 
noted that E x [p* , pi] is an implicit functional of the spin 
densities p a {v). Of course, correlation can be considered 
as well in a similar framework. In this kind of functionals 
the calculation of the exchange-correlation potential re- 
quires the OEP method . 22 ' 23 For a review of the method 
the reader is referred to Refs. [28ll29l[30l . Restricting our- 
selves to the EXX case, the OEP method leads to an 
integral equation for the exchange potential, which can 
be written in compact form as 

^fe(r)^(r)+c.c.)=0, (11) 
i=i 

where the so-called orbital shifts are defined as 
V* (7 (r) = J d 2 r'^v')[v x<y {v')-u X3a {v')]G%{v',v). 

(12) 



Here G s ja is the Gr een function of the KS system, 

^(r',r)= ± *k(r> to (r) , (13) 

k — l J 

and 

U x ja(Y) = — . (14) 

In a series of steps ; 30 ' 31 the OEP equation as given in 
Eq. (fTT|) can be transformed to 
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r(r) = 



2p.(r) 



E 



\fja{r)\ 2 (u xja (r) + (v xja - u xja )) - V • (V>* CT (r)V^ CT (r)) 



(15) 



where 



and 



= / d 2 r V * a {v)v xa {v) V3a {v) , (16) 



d 2 r ip* j<7 {r)u xja {r)<p ja {r) . (17) 



The OEP equations can be solved iteratively, simulta- 
neously with the corresponding KS equations, in a self- 
consistent fashion. 

Due to the presence of the unoccupied KS orbitals in 
the definition of the orbital shifts [see Eqs. (fT2"j) and [fT5])]. 
the full numerical solution of the OEP integral equation is 
nontrivial. Of course, one may take advantage by spec- 
ifying it for a particular kind of systems. In the origi- 
nal paper, 23 solutions were presented for atomic systems 
with spherical symmetry. Much later, the OEP equation 
has also been solved for systems with lower symmetry 
such as molecules ; 32 ' 33 solids metallic surfaces^ and 
quasi two-dimensional electron gases at the interface of 
two different semiconductors^ In addition, an iterative 
algorithm for the solution of the OEP equation based on 
the orbital shifts has been implemented . 36 ! 37 Recently, 
progress has been made in studying physically relevant 
examples of non-collinear magnetism ) 38 ' 39 and open-shell 
systems in the relativistic limiti 40 ' 41 On the other hand, 
one may approximate the full OEP in order to save com- 
putational effort. Along this line, the KLI approach 2 ^ has 
turned out to be rather accurate in many situations. In 
the KLI approximation, the terms containing the orbital 
shifts on the r.h.s. of Eq. (jTSJ) are neglected completely. 30 

For the analysis which follows in the next sections, we 
may rewrite Eq. (fT5|) as 



V x <r( v ) 



^ (r)+A ^ LI (r) + A^ (T s (r), 
^ L (r) + Az,° EP (r), 



(18) 



where 



= T-7T E l^>(r)|V;a(r) + C.C. (19) 
^ ' i=i 



is the so-called Slater potential; 

N, 



+ c.c. 



is the correction to the Slater potential added by the KLI 
approximation of the OEP solution; and 



-1 



Av os (r) = — V 



^ CT (r)| 2 V-(^ a (r)V^ CT (r))- 



-c.c. 



(21) 

is the correction to the KLI solution, which adding the 
contribution from the orbital shifts, provide the full OEP 
solution. Thus, the definition of Aw° EP (r) in Eq. HU) 
becomes clear. We remind that the asymptotic limit of 
v X (t(y) in Eq. (fT8|) . for finite systems, is given by the first 
term, v^(r — > oo) — > — 1/ n 14 ' 30 

We point out that the derivations above are fully inde- 
pendent of the dimensionality. In orbital functionals the 
dimensional character is ensured from the structure of the 
KS orbitals themselves, and no further assumptions are 
made. Therefore, the OEP (and KLI) schemes presented 
above apply equally well to 3D and 2D. Below, how- 
ever, we focus on a recent approximation to Av x ® F pro- 
posed by Becke and Johnson, 25 (BJ) which is dimension- 
dependent by construction. In 3D, the BJ approximation 
has received considerable attention. 42 - 43 - 44 ' 45 - 46,47 - 48 ' 49 



III. BECKE-JOHNSON APPROXIMATION 

A. Three dimensions 

For 3D systems, Becke and Johnson have proposed to 
approximate Ai)°^ p (r) by Au^(r), defined as^ 5 - 



1 
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[r.(r)l 







with 



r CT (r) = ^|V^(r)| 2 



(22) 



(23) 



k=l 



(20) 



being (twice) the spin-dependent kinetic-energy density. 
This approximation was found by seeking a simple ex- 
pression having the following properties: (i) it is (possi- 
bly) invariant with respect to unitary orbital transforma- 
tion; (ii) it provides an exact treatment of any (ground 
state) hydrogenic atom; (iii) it has a step-like structure 
characteristic of Av x ^ p (r) in multi-shell atoms; 28 ' 29 ' 30 
(iv) it is exact in the 3D uniform electron gas limit, 
Av° CT EP = Av^ = [3p CT /(47r)] 1/3 , with constant spin den- 
sities. The BJ potential as defined in Eq. (|22|) has been a 
good approximation for several atoms^ and it has given 
valid estimations for semiconductor band gaps i ' 
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B. Extension to two dimensions 

In the framework of the BJ approach, we may suggest 
the following approximation to be valid in 2D, 



A« Xff (r) = 



2^/2 
~3tP 



P<x( r ) 



(24) 



where the coefficient 2y/2/(3n) is chosen to satisfy the 
exact 2D uniform electron gas limit 



Av 



OEP 



= Av xrr = 



,V2 



(25) 



with constant densities. However, as shown analytically 
and also numerically in the next section, this form serves 
only as the starting point. We need to introduce two 
completely general modifications, based on an additional 
set of physical requirements. 

First, we observe that property (i) given in the pre- 
vious section is satisfied only under the restriction of 
real- valued orbitals. In fact, it is easy to see that under 
gauge transformation for both real and complex-valued 
orbitals the BJ contribution to the exchange potential 
is not gauge invariant. The relevance of the gauge in- 
variance requirement for the meta-generalized-gradient 
approximation has been analyzed in some detail by Tao 
and Perdew^S, We expect that the same criterion should 
apply to any approximation for v xa (r) which depend ex- 
plicitly on T CT (r). We can ensure the gauge invariance by 
replacing r CT (r) wit h 19 ' 50 



(26) 



where 



N„ 



jp^(r) = \ % £ bL(r) (VptoCr)) - (V^LW) Wtar(r)] 

k=l 

(27) 

is the spin-dependent paramagnetic current density. 

Second, property (ii) given in the previous section is 
not valid for all one-electron systems. As we show be- 
low, this also relates to the wrong (divergent) behavior 
of the approximation in the asymptotic region. As it has 
been shown in 3D, an ad- hoc solution may be found*^ for 
specific cases, but still it would be preferable to have a 
general cure. Hence, to obtain an expression which van- 
ishes in the case of a generic single-particle state, we add 
an additional term 



7v(r) — > D a {r) = f CT (r) - 



In conclusion, in 2D we obtain 



1 (MIT 

4 Pa(r) 



Av xa (v) 



(r) 



2V2 

^T 1 



Ar(r) 



(28) 



(29) 



The present approximation for the total corrected 2D ex- 
change potential reads then 



tw(r) = v^(r)+Av] 



SL 



(r) 



2V2 \D a (r) 



3tt 



with 



4 pa{r) 



P<y{r) 



Pa{r) 



(30) 



(31) 



Eq. (|30|) is the key expression of the present work, and in 
the following we will discuss its properties, first analyt- 
ically and then numerically. First, it is easy to see that 
D a (r) = for any one-electron system (per spin channel) , 
as it should be [property (i)]. To see this, it is enough to 
note that for this particular case p a {r) — |</3 CT (r)| 2 , with 
¥<j{ r ) = V Pa(r)e ie ( r \ with no loss of generality. Sub- 
stituting this in Eq. (j3Tj) . the result Do-(y) = follows 
at once, for arbitrary functions p a (r) and 6(r). Second, 
in the asymptotic limit, D a (r — > 00) — > 0, for any fi- 
nite system. To see this, it is enough to note that for 
any finite system, all ingredients of D a (r) in Eq. (|3"Tj) 
become dominated by the contribution of the highest 
occupied orbital. As a consequence, in the asymptotic 
region the system becomes effectively one-electron like, 
and the result D a (r — > 00) — > applies. This guar- 
antees that v X(T (r) has the correct asymptotic behavior 
given by i£^(r). Besides correcting the potential in the 
asymptotic region, we also expect improvements at finite 
r. We find this reasonable due to the fact that D a (r) 
is the main ingredient of the so-called electron localiza- 
tion function ) 56 i 57 ' 58 which is able to deal with situations 
where the multi-shell structure of atoms becomes rele- 
vant, as well as in the region of interest for chemical 
bonding. 

Then, it is reassuring to note that one can justify the 
choice of the power a in 



A^(r) 



2^2 
~3V 



Ax(r) 



(32) 



by observing that the resulting potential scales linearly 
only for a = 1/2, similarly to the exact exchange 
potential^ Under uniform scaling of the coordinates, 
r — > Ar (with < A < 00), and for the norm-preserving 
many-body wave function, the 2D density scales quadrat- 
ically with A, p^( r ) ^ 2 Pa(^ r )- This leads to the 
result that the KS orbitals in 2D are seen to scale 



A <^ feCT (Ar). Thus, 



Vp^(r) -> A 3 V Ar p CT (Ar), and j^(r) 
stituting these relations to Eq. ([32 
relation 



r» - A 4 T(j ( Ar) , 
-> \%a{\r). Sub- 
yields the scaling 



AvZ\r) -» X 2a AvZ(Xv) 



(33) 



which fulfills the exact linear scaling constraint only if 
a = 1/2. 
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Finally, we remind that the calculation of Av xa (r) is 
computationally cheap in comparison with the full OEP 
method or the KLI approximation, which is naturally the 
case also for the conventional A-y^(r). The calculation 
of the Slater part of the total potential is actually still 
costly. For this part, we have previously provided ac- 
curate and numerically simple approximations in 2D by 
considering the properties of the exchange holei 13 i 14 



IV. RESULTS 

A. Single-electron states of a harmonic oscillator 

First, we consider the single-electron non- interacting 
eigenstates in a 2D isotropic harmonic oscillator. The 
external confining potential is given by VQ a {r) = uj 2 r 2 /2, 
where uj is the oscillator strength. The use of this sys- 
tem is motivated by the frequent use of such a confining 
potential when modeling 2D quantum dots^ (see also be- 
low). The eigenstates, with the associated densities, are 
given by 



Mr,0) = U{r) e- r2 ' 2 
Pni(r) = \<fini(r,9) 



fni(r)e- r , 



(34) 
(35) 



respectively. Here r and 9 are the usual radial and 
angular 2D polar coordinates, n = 0,1,2,... and I — 
0, ±1,±2, ... are the radial and angular quantum num- 
bers, and /ni( r ) is a radial function, related to Laguerre 
polynomials. For simplicity, we do not consider the spin 
index in this analysis. The corresponding contribution 
to the kinetic energy density is given by 



Tni(r) = |W„,(r,0)| 5 



dfni jr) 
dr 

l 2 Pni(r) 



rfni{r) 



Following Eq. (|24|) . we obtain 



df n i(r)/dr 
fm(r) 



Pm(r) 



(36) 



(37) 



Clearly, this contribution to Av xa (r) diverges linearly for 
r — > oo, after summing over all occupied states (all show- 
ing the same divergent behavior), and after taking the 
square root of the sum. On the other hand, this divergent 
behavior is absent in our corrected exchange potential in 
Eq. ([3H]) . as it cancels exactly with the second term of 
Ar(r) in Eq. QSJ. 



B. Many-electron quantum dots 

Next we focus on the same system as in the previ- 
ous section but consider N interacting electrons. Within 
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FIG. 1: (Color online) Exchange potentials for a two-electron 
quantum dot (singlet). The results are calculated from the 
analytic density. Av X a and Av^ a are defined in Eqs. (J24J) and 
(|29l) . respectively. 



the effective-mass approximation for electrons in semi- 
conducting host material such as GaAs, this is the most 
common model for 2D quantum dotsi^ 

First, we consider the singlet solution of N = 2 
and u> — 1 a.u., which is known analytically^ In 
this case, the exact exchange potential is equal to the 
Slater potential (and minus half of the Hartree poten- 
tial). Therefore, any correction to the Slater poten- 
tial should be zero. For the kinetic- energy density we 
find T a (r) = (V/9 (T (r)) 2 /(4p cr (r)), and the non-corrected 
Becke- Johnson-type approximation in Eq. ([24]) becomes 
Av xa (r) — ^/2Vp a (r)/(37rp CT (r)), where the spin density 
has an analytic expression 



1 + r 2 /2) + -V7T< 



-3r 2 /2 



7T(V27T + 3) 

I (r 2 /2) + r 2 h{r 2 /2) + r 2 h(r 2 /2)] }, (38) 



deduced from the analytic wave function^ 2 - Here Iq and 
I\ are the zeroth and first-order modified Bessel func- 
tions of the first kind, respectively. Figure [T] shows that 
Av xa (r) (dotted line) has clear divergent behavior, mak- 
ing the total exchange potential (dashed line) largely dif- 
ferent from the exact result (thick solid line), even in the 
small-r regime. It should be noted that also the LDA ex- 
change potential considerably deviates from the exact re- 
sult. On the other hand, our corrected scheme in Eq. (|3"0|) 
yields the exact result: Av xa (r) = as discussed above. 

The divergent behavior of Av xa (r) may be easily ex- 
tracted from the interacting density given in Eq. (|38[) . 
For this, one should use the asymptotic expansion of the 
modified Bessel functions: l v (z — > oo) — * e z / \J2t\z + ... 
(Ref. |H). Using this in Eq. IJ38J leads to p a (r) -> 

r 2 e- r2 /7r(3 + v / 2^ : ), and Vp a (r) -> -2r 3 e - r2 /Tt(3+V2w) 
in the asymptotic region (r — > oo). Replacing everything, 
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FIG. 2: (Color online) Exchange potentials for a fully spin- 
polarized three-electron quantum dot in external magnetic 
field of B = 5 T. Terms with Av^ a and Av^ a refer to the 
corrected approximation [Eq. (|30|l ]. with and without the ex- 
plicit current term in Eq. pip , respectively. The scales for the 
vertical axis referring to the potentials and density (negative 
and positive values, respectively) are different. 




2 4 6 8 

r (eff.a.u.) 



FIG. 3: (Color online) Same as Fig. [2]but for ten electrons at 
B = 8 T. 

we obtain that Av xa (r — » oo) — > 2v / 2r/(37r) ~ 0.300 r 
in agreement with Fig. [TJ It is interesting to note that 
the leading term to the divergent contribution obtained 
for the interacting system agrees exactly with the non- 
interacting estimate given above in Eq. (|37[) . The reason 
is that the density is rather small in the asymptotic re- 
gion, and then the modifications introduced by interac- 
tions to the bare 2D harmonic potential are also rather 
small. As a consequence, the interacting densities ap- 
proach the non-interacting limit asymptotically. 

In Figs. [5] and [3] we show the exchange potentials 
v xa for quantum dots containing three and ten interact- 
ing electrons, respectively. In both cases the oscillator 



strength is ui = 0.4217 in effective atomic units (eff.a.u.) 
with the material parameters for GaAs: m* = 0.067 m e 
and e = 12.4 eo- Here we have performed self-consistent 
KLI calculations by applying the octopus DFT code, 54 
and use the resulting KS orbitals as inputs in the ap- 
proximations given above. Both dots are exposed to 
an external, perpendicular magnetic field with a suffi- 
cient strength to achieve full spin polarization (S — N/2) 
and occupation to eigenstates of consecutive angular mo- 
menta from I = to I = — N + 1. The resulting 
"maximum-density droplet"— has a smooth electron den- 
sity (see the upper panels of Figs. [5] and [3]). The cor- 
responding exchange potentials calculated with different 
approximations have several interesting features listed 
below. 

• Similarly to the previous examples, the straightfor- 
ward BJ extension, i.e., Slater potential combined 
with Av xa in Eq. (|24p . shows divergent behavior. 

• Our approximation Av^fr) in Eq. (f3"0"]) leads to 
good overall agreement with the KLI potential, es- 
pecially at large r. The agreement improves as a 
function of N. 

• The LDA exchange potential is overall worse than 
our approximation - not only in the asymptotic 
region showing the obvious exponential decay - but 
also close to the center of the system. 

• The current-dependent term introduced in 
Eq. (|26p . which enforces the gauge invariance, is 
crucial in finding the correct behavior. Without 
this term (Av^, a ) the performance is very poor at 
high magnetic fields (see in particular Fig. [3] with 
B = 8T). 

• Close to the center of the maximum-density 
droplet, which locally resembles the 2DEG with the 
electron localization function « 1/2 (Ref. |58[). the 
present approximation satisfies with a good pre- 
cision the well-known exact relation for the ho- 
mogeneous 2DEG between the exchange potential 
and the Slater potential: v xa = (3/4)wfJ; (see the 
crosses in Figs. [3] and [3|). For the calculation of the 
cross positions, the numerical value of v^(r = 0) 
has been multiplied by 3/4. Note that this rela- 
tion is not well satisfied between v x ^ A (r = 0) and 
v xa( r = 0): because the quantum-dot system is fi- 
nite, and globally not homogeneous. 

Figure [1] shows the exchange potentials for a closed- 
shell quantum dot with 20 electrons. Again, the straight- 
forward BJ extension to 2D diverges, whereas the cured 
approximation qualitatively agrees with the KLI solu- 
tion. It should be noted that the shell structure is better 
reproduced by the present approximation than by the 
LDA. In this system, the (3/4)-relation is not well sat- 
isfied (see the cross in Fig. @|, since the system proper- 
ties in the core region are far from those of the homo- 
geneous 2DEG, i.e., the electron localization function is 
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FIG. 4: (Color online) Exchange potentials for a closed-shell 
quantum dot with 20 electrons. 

varying. 58 Also note that for this closed-shell system, the 
spin-dependent paramagnetic current density vanishes, 
and then Av° a (r) = Av° a (r). 

Finally, we point out that the present scheme, in par- 
ticular the modifications introduced in Eqs. (j26|) and 
(|28p. can be considered and reworked also in 3D. This 
has been shown to lead to an approximation for the 3D 
exchange potential that gives good results in systems 
possessing orbital currents (e.g., open-shell systems), in 
the presence of external magnetic or electric fields, or in 
"non-atomic" systems such as the Hooke's atom, i.e., the 



3D counterpart of the harmonic quantum doti^ 



V. CONCLUSIONS 

To achieve an accurate and computationally conve- 
nient approximation to the exchange potential in two 
dimensions, we have applied and extended the formal- 
ism of Becke and Johnson commonly used in three di- 
mensions. We have found that direct extension of the 
formalism leads to divergent behavior of the exchange 
potential in the most standard two-dimensional appli- 
cations. We have corrected the scheme by introducing 
a current-dependent term ensuring the gauge invariance 
and a density-gradient term capturing the disappearance 
of the correction to the Slater contribution for a generic 
single-particle state. The resulting approximation for the 
exchange potential agrees well with the Krieger-Li-Iafrate 
potential in various quantum dots up to high magnetic 
fields and thus to high current densities. The correct 
asymptotic behavior is recovered and the shell-structure 
is well described. Similar strategy is expected to lead 
to good results also in three-dimensional applications in 
external magnetic and/or electric fields. 
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